Formation of quark phases in protoneutron stars: the transition from the 2SC to the 

normal quark phase 
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We study the process of formation of quark phases in protoneutron stars. After calculating the 
phase transition between nucleonic matter and the 2SC phase at fixed entropy and lepton fraction, 
we show that an unpairing transition between the 2SC phase and the normal quark phase occurs for 
low lepton fractions. We then calculate the process of diffusion of neutrinos in protoneutron stars 
and show that for intermediate values of the mass of the star, the deleptonization triggers the phase 
transition between the two quark phases after a temporal delay of a few seconds. In less massive 
stars instead only the normal quark phase is formed at the end of the deleptonization stage. We 
also discuss the possible astrophysical implications of our scenario. 



I. INTRODUCTION 

In the last years, a new very exciting theoretical discov- 
ery in QCD concerns the low temperature and high den- 
sity region of the phase diagram of strongly interacting 
matter where quarks form a state of color superconduc- 
tivity. In particular, the so called Color-Flavor-Locking 
(CFL) phase in which all quarks are paired, was shown 
to be the ground state of strongly interacting matter at 
asymptotically high densities [1]. In turn, this opens the 
possibility that this new quark phase appears in the core 
of neutron stars and affects their properties. A number 
of studies has been done to investigate the structure of 
the QCD phase diagram for the conditions that are re- 
alized in neutron stars: the requirements of charge neu- 
trality and equilibrium with respect to weak interactions, 
which must be fulfilled in compact star matter, and the 
finite value of the mass of the strange quark split the 
values of the chemical potentials of the different quark 
flavors thus allowing for the existence of many color su- 
perconducting phases in the phase diagram as the 2SC 
phase, the normal (unpaired) phase, the gapless phases 
and the crystalline phase [2j-|6j|. An important question 
concerns when and how these quarkphases are eventually 
formed in a neutron star. In 0, 043 it was argued that 
while the 2SC phase or the normal quark phase might be 
formed already during the protoneutron star stage i.e. 
when matter is hot and lepton rich, the CFL phase could 
appear only later when the star has almost completed its 
deleptonization through the emission of neutrinos. The 
possibility of forming different quark phases during the 
protoneutron stars evolution could have very intriguing 
observable consequences as proposed in (lol - [T3| in con- 
nection with gamma-ray-bursts. 

To date, there are no detailed numerical simulations of 
the process of formation of color superconducting phases 
in protoneutron stars: such a calculation is indeed very 
complicated because, besides the equation of state, it in- 
volves also the computation of the neutrino mean free 
path in the different quark phases and it requires to solve 
numerically the transport equations describing the diffu- 
sion of neutrinos within the protoneutron star matter. 



A first preliminary study was presented in [l4|, where, 
besides the calculation of the neutrino mean free path 
in color superconducting quark matter, a simple mod- 
eling of the neutrino cooling was presented in order to 
discuss the possible observable signatures of the forma- 
tion of a color superconducting phase in the star. The 
scenario consists of a second order phase transition from 
the unpaired phase to a superconducting phase: during 
the cooling, the temperature decreases and at some point 
reaches the critical temperature and the peak of the spe- 
cific heat associated with the phase transition is met. The 
authors argue that the phase transition would lead to a 
slower cooling of the star for temperatures close to the 
critical temperature with possible effects on the neutrino 
luminosity. This possibility could be indeed realized in 
protoneutron stars if the color superconducting gap is 
of the order of a few tens of MeV which lead to criti- 
cal temperatures comparable with the temperatures of 
protoneutron stars. 

Here we will investigate a different scenario: at birth 
a protoneutron star has a lepton fraction of the order of 
Yl = 0.4 with a corresponding electron/proton fraction 
of Y e = Y p ~ 0.3 at the center of the star. The sub- 
sequent deleptonization causes a decrease of lepton and 
also electron fractions: consequently also Y p decreases 
and matter becomes more and more isospin asymmetric. 
If quark matter is present in the core of the star at birth, 
it is likely to be in the 2SC phase because matter is not 
very isospin asymmetric i.e. the number of up and down 
quarks are quite similar and the 2SC pairing is favored. 
As the deleptonization proceeds, the gradual increase of 
the isospin asymmetry causes a stress on the Cooper pairs 
and a first order phase transition from the 2SC to the un- 
paired quark phase is trig gered at some critical value of 
the asymmetry 0, Il5l4l8j | Remarkably, the phase tran- 
sition between the superconducting state and the nor- 
mal state in asymmetric fermionic systems has attracted 
much attention in the last years also in connection with 
cold atoms experiments (see for instance [l9|, Ho[). By 
solving the neutrino diffusion equation coupled with a 
quasi-static temporal evolution of the structure of the 
star, we will show here that such a phase transition be- 
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FIG. 1: Pressure as a function of the baryon density for nu- 
cleonic matter for values of Yl ranging from 0.4 (thick line) 
to 0.12 (dashed line) with intervals of 0.04 (solid lines). At 
lower lepton fractions, which implies also lower electron and 
proton fractions, the nucleonic equation of state is stiffer be- 
cause of the increasing contribution of the symmetry energy 
to the EoS. 



tween quark phases could indeed occur in protoneutron 
star matter and possibly gives some observational effects 
in the neutrino signal. 

The paper is organized as follows: in Sec. II we cal- 
culate the Equation of State (EoS) for protoneutron star 
matter by including the pure quark and nucleonic phases 
and the mixed phase. In Sec. Ill we present and discuss 
the results of our numerical simulation of the neutrino 
diffusion and finally in Sec. IV we draw our conclusions. 



II. EQUATION OF STATE OF 
PROTONEUTRON STAR MATTER 

Matter in protoneutron stars is hot and lepton rich. 
At birth, the initial conditions correspond, to good ap- 
proximation, to uniform values of the entropy per baryon 
S/N = 1 and the lepton fraction Yl = 0.4 [2l|,[22j]. Since 
the deleptonization reduces the value of Yl , we will cal- 
culate the EoS for Yl varying within the range 0.04 — 0.4. 
As we will explain in the next section, we disregard here 
the effect of reheating of the star during deleptonization 
and we will therefore keep S/N fixed during the simu- 
lation. To describe the phase transition from Nucleonic 
Matter (NM) to Normal Quark matter (NQ) or 2SC mat- 
ter, we use, as customary, two models: for the low den- 
sity nucleonic phase we adopt the relativistic mean field 
model with the parameterization TM1 [23[ and for the 
high density quark phase a modified bag model which in- 
cludes also color superconductivity in the 2SC phase used 
in [l5j]. We limit our discussion here to two flavor quark 
matter, with massless up and down quarks and we fix the 
model parameters as follows: the intermediate value of 
the diquark coupling is considered Go — 3 /AGs where 
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FIG. 2: Pressure as a function of the baryon density for the 
nucleonic matter - quark matter phase transition for values 
of Yl ranging from 0.4 (thick line) to 0.12 (dashed line) with 
intervals of 0.04 (solid lines). The dots on the dashed line 
indicate the onset and the end of the mixed phases of nucle- 
onic matter and normal quark matter and nucleonic matter 
and 2SC phase. At large values of Yl , the quark phase corre- 
sponds to the 2SC phase. At lower Yl, the mismatch between 
the chemical potentials of up and down quarks is such that 
the pairing is broken and the normal quark phase is obtained. 



Gs is the scalar coupling and the bag constant B is fixed 
to 155 MeV 1 / 4 . We also include, in an effective manner, 
the corrections to the quark pressure due to perturba- 
tive strong interactions [24], [25|: the constant c, which 
simulates perturbative interactions, is fixed to a value 
of 0.05 for which we obtain a maximum mass for cold 
and catalyzed hybrid stars of 1.94M , in agreement with 
the very recent measurement of a compact star having a 
mass of 1.97 ± 0.04M© [26[. For this choice of param- 
eters, the onset of the phase transition occurs at twice 
nuclear saturation density for the initial lepton fraction 
and entropy configurations. The phase transition is com- 
puted by using the Gibbs construction as described in 
[2ll I22I [27rl29l | with the conditions of global conservation 
of baryonic number, lepton number and electric charge. 
The resulting extended mixed phase is obtained by solv- 
ing the following equations: 



(1 - X )nc +xn2-n e = 
n e + n v = Y L n B 
(1 - X)s H + Xs Q = SlNn B 
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(1) 

(2) 
(3) 
(4) 
(5) 
(6) 



where p\ [i B , [i l c , T 1 , n' c , s l are the pressure, baryon 
chemical potential, charge chemical potential, lepton 
chemical potential, temperature, electric charge density 
and entropy density of the hadronic (i — H) and the 
quark phases (i = Q); n e and n v are the electron and 
electron neutrino densities, ub is the baryon density and 



X is the volume fraction of the quark phase. It has been 
recently pointed out that another possible hadron-quark 
mixed phase could appear in protoneutron stars if lo- 
cal charge neutrality and global lepton number conserva- 
tion are imposed [30| ; we will not consider this possibility 
here. 

In Fig. [1] we show the EoS (pressure vs baryon den- 
sity) of nucleonic matter for different values of Yl . Notice 
that by decreasing Yl, the EoS becomes stiffen despite 
the fact that at lower Yl the contribution to the total 
pressure of leptons is smaller, the lower values of Y e im- 
plies a larger isospin asymmetry of nucleonic matter. The 
effect of the symmetry energy thus dominates over the 
reduced leptonic pressure and the EoS is stiffer. This a 
well known fact and it also implies that the maximum 
baryonic mass of protoneutron stars is smaller than the 
maximum baryonic mass of cold neutron star (if nucle- 
onic matter is considered) [3lj . 

Let us discuss now our results for the mixed phase. 
We show the EoS in Fig. 2: at a large lepton fraction, 
Yl = 0.4, after a low density nucleonic phase, the 2SC- 
NM mixed phase starts and at high density the pure 2SC 
phase takes place. By decreasing the value of Yl, at 
some point, close to the onset of the mixed phase, the 
2SC pairing cannot take place anymore because of the 
large mismatch between the up and down chemical po- 
tentials (the isospin density increases as Yl decreases) 
and the normal quark phase starts to form. Notice that 
the transition from the 2SC phase to the normal phase 
occurs close to the onset of the mixed phase and moves to 
higher densities as Yl is further decreased. This is due to 
the fact that the local isospin density of the quark phase 
is large at the onset of the mixed phase and gradually 
decreases as the volume fraction increases, as shown in 
Refs. [TH, HI, HII- The largest stress on the up and down 
chemical potentials is thus realized close to the onset of 
the mixed phase: the hadron quark phase transition in- 
duces also a phase transition between quark phases. To 
calculate the transition between the NQ and the 2SC 
phases within the mixed phase we use the following sim- 
ple prescription: we start with the high density pure 2SC 
phase, x — lj an d we decrease the density until we reach 
the border of the mixed phase with the NM obtained by 
solving the system of Eqs. 1-6; we solve the same system 
at decreasing values of the baryon density and x and ob- 
tain the solutions fi* B , [i* c , \x* L and T*; then we accept 
the solution of the system only if the pressure of the nor- 
mal quark phase computed at the same [i* Bl [i* c , fi* L and 
T* is lower than the pressure of the 2SC phase. When 
this condition is not fulfilled anymore we choose the NQ 
phase as the quark component of the mixed phase. This 
is a simple way to describe the unpairing phase tran- 
sition between the 2SC and the NQ phases (similar to 
a Maxwell construction) within the mixed phase with 
NM. Actually, since this phase transition is also of first 
order one should introduce two volume fractions associ- 
ated with the two quark phases xnq and X2SC and solve 
a more complicated system of equations which provides a 
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FIG. 3: Temporal variation of the neutrino density as a func- 
tion of the enclosed baryonic number for a neutron star of to- 
tal baryonic number Nb = 2.29 x 10 57 . The initial central neu- 
trino density at the center of the star, n„ = 1.88 x 10 -2 fm -3 . 
The thick line corresponds to the initial configuration, t = 0, 
and the dashed line to the configuration at t — 10s. The 
thin solid lines stand for the intermediate configurations with 
temporal steps of Is. 



mixed phase made of three components. A similar calcu- 
lation has been done, for cold and beta stable matter, for 
mixed phases composed by quark color superconducting 
phases [34j |. but calculations for protoneutron star mat- 
ter, also including the nucleonic phase, have never been 
performed. Our procedure, while simple, is sensible since 
the NQ-2SC phase transition is not a strong first order, 
i.e. the jump of the density is not large when going from 
one phase to the other one (the density of the 2SC phase 
is larger than the one of the NQ by an additive term 
scaling as (A//i) 2 , with A ~ 100 MeV being the super- 
conducting gap) . Therefore the three components mixed 
phase would actually cover a small interval of density. 
We consider a full calculation with the three components 
as an interesting future development of this work. 



III. NEUTRINO DIFFUSION IN 
PROTONEUTRON STARS 

The evolution of a protoneutron star can be schemat- 
ically divided into two separate stages [31| : the delep- 
tonization, during which the trapped neutrinos diffuse 
and are gradually released from the star, and the cool- 
ing, during which the entropy of the star reduces by the 
emission of free streaming neutrinos. During deleptoniza- 
tion, the star is actually reheated by the interactions of 
neutrino with baryonic matter and the initial S/N ~ 1 
raises to a value of ~ 2. This is for us the most inter- 
esting stage because Y e and Y p also decrease during it 
thus leading to a sizable change of the chemical compo- 
sition of the star which, as we will show, is responsible 
for the formation of different quark phases. The stan- 
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FIG. 4: Temporal variation of the baryon density as a func- 
tion of the enclosed baryonic number for a neutron star of 
total baryonic number Nb = 2.29 x 10 57 . The thick line cor- 
responds to the initial configuration t = and the dashed 
line to the configuration at t = 10s. The thin solid lines stand 
for the intermediate configurations with temporal steps of Is. 
A small decrease of the central density occurs during delep- 
tonization. 




FIG. 5: Temporal variation of the pressure as a function of 
the enclosed baryonic number for a neutron star of total bary- 
onic number Nb = 2.29 x 10 57 . The thick line corresponds 
to the initial configuration t = and the dashed line to the 
configuration at t = 10s. The thin solid lines stand for the in- 
termediate configurations with temporal steps of Is. A small 
decrease of the central pressure occurs during deleptonization. 



dard technique to simulate the evolution of protoneutron 
stars consists of solving two partial differential equations 
associated with the transport of lepton number and the 
transport of energy within the so called diffusion approx- 
imation (35M38| . These equations are coupled with the 
Tolman-Oppenheimer-VolkofF (TOV) equations to take 
into account the mechanical readjustment of the star dur- 
ing the evolution. Indeed, the time scales of diffusion are 
much larger than the dynamical time scale of the compact 
objects and therefore the evolution of the structure of the 
star follows quasi-static equilibrium configurations. In 
the cases in which the deleptonization drives a collapse of 
the star, because of the formation of exotic phases for in- 
stance, the hydrodynamical evolution must be also taken 
into account as done in [3^] or more recently in [4l[ 
where for the first time consistent simulations from the 
core collapse of the supernova to the cooling epoch of the 
newly formed neutron star have been performed. 

The crucial ingredients for such calculations are the 
EoS and the cross sections of all the possible interac- 
tions of neutrinos with the matter of the star. Different 
studies have addressed the possibility to form "exotic" 
phases during the evolution of a protoneutron star as 
hyperons [37], HH, kaon condensates [42[ or finally nor- 
mal quark matter (43j . A similar calculation taking into 
account color superconducting phases is still missing and 
this is the problem that we want to address in this paper. 
Since we want here to provide a qualitative description 
on how the unpairing transition could occur in a compact 
star, we will adopt the following simplifying assumptions 
in the treatment of the neutrino transport in protoneu- 
tron stars: we will consider only the first seconds of the 
stage of deleptonization when neutrinos are still degener- 



ate and we will assume the entropy profile to be constant 
during the evolution (S/N = 1). These assumptions lead 
to several simplifications: we will not solve the equation 
for the energy transport but just the one describing the 
lepton number transport coupled with the TOV equa- 
tions; when neutrinos are degenerate the dominant con- 
tribution to the mean free path is given by the absorption 
processes of neutrinos by non-degenerate neutrons in NM 
and by degenerate quarks in the quark phases [U HiJ ■ 
Moreover, also regarding the EoS to be used as input for 
the simulations we will use a two-dimensional table with 
the baryon density and the lepton fraction as indepen- 
dent variables and the ratio S/N will be kept constant. 
Because of our simplifying assumptions, our calculation 
has two important limitations: the estimate for the tim- 
ing of the phase transition must be taken only as an order 
of magnitude estimate; indeed, because of the reheating 
that we are neglecting here, all the cross sections of neu- 
trinos would be enhanced by the increase of the temper- 
ature and therefore we expect the global evolution of the 
system to be slower than what we find here. Secondly, 
we cannot provide here an estimate of the neutrino lumi- 
nosity which is obtained by solving the energy transport 
equation. 

Let us discuss now the equations that we solve to sim- 
ulate the deleptonization [371 ]: 
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FIG. 6: Temporal variation of the neutrino density as a func- 
tion of the enclosed baryonic number for a hybrid star of bary- 
onic number Nb = 2.29 x 10 57 . The initial central neutrino 
density, n„ = 3.44 x 10 -2 fm -3 . The thick line corresponds 
to the initial configuration t = and the dashed line to the 
configuration at t = 10s. The thin solid lines stand for the 
intermediate configurations with temporal steps of Is. 



da 
dr 
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dr 



yl — 2m /r 
m + 4irr 3 p 
r 2 — 2mr 



(10) 
(11) 



where r is the radial distance in spherical coordinates, 
r = y 7 ! — 2m/r, and are the general relativity correc- 
tions to the diffusion equation, A is the spectral average 
of the electron neutrino mean free path, p, e, m and a 
are the pressure, the energy density, the enclosed gravi- 
tational mass and the enclosed baryonic number (we use 
the gravitational units G = 1 and c = 1). The first 
equation describes the diffusion of neutrinos driven by 
the neutrino density gradient (i.e. within the diffusion 
approximation ) and the last four equations describe the 
mechanical equilibrium of the star, the TOV equations. 
The initial condition and the boundary conditions are 
specified as follows: for the initial condition we use an 
EoS with constant S/N = 1 and Yj, = 0.4; on the bound- 
aries, r = and r = R (where R is the radius of the star 
which, numerically, corresponds to a very small value of 
the pressure), the neutrino flux F(r) — — ^-fer- is taken 
as F(0) — and F(R) — fn Vc where / is a constant 
varied between 0.1 — 0.6. The boundary condition at 
r = R basically imposes that the flux of neutrinos is 
a certain fraction / of the flux that one would have if 
neutrinos were free streaming [381 ] . We also need ini- 
tial conditions for the TOV equations: the central pres- 
sure is chosen arbitrarily, p{r = 0) = p c (different val- 
ues of p c lead to different values for the total baryon 
number of the star), m(r = 0) = 0, a(r = 0) = and 
4>(r = R) = 7}log(l — 2m/ V). Finally, concerning the 
mean free paths, during the initial stage of the evolu- 
tion, the absorption processes are the dominant ones and 



thus we take for NM the mean free path associated with 
the absorption of degenerate neutrinos by non-degenerate 
neutrons 1311: 



n n a (l + 3g A ) 



E v 



(12) 



where n n is the neutron number density, cro = 1.76 x 
10 -44 cm 2 , gA = 1.257, m e is the electron mass and 
E v , the neutrino energy. This formula must be actually 
corrected by a factor of 3 — 10 to take into account the 
degeneracy and Fermi liquid corrections [3lj ] , we will use a 
factor of ten in this paper. Similarly for quark matter we 
consider the absorption processes of degenerate neutrinos 
by degenerate down quarks 
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2G>3( 10/i 2 + 5 Mu/Xe + M 2)((£„ - M „)2 + ^2 T 2) 

This formula holds for normal quark matter and it can be 
used also for the blue quarks of the 2SC phase which are 
indeed unpaired. We therefore multiply X q by a factor 
of three within the 2SC phase to take into account that 
the absorption process can take place only if blue quarks 
are involved. For the paired red and green quarks of the 
2SC phase, the corresponding mean free path is larger 
[l4j then the one associated with the blue free quarks 
and thus we can safely neglect them in our calculation. 
Finally, since we study the regime in which neutrinos 
are still degenerate, the spectral average of the mean 
free path is dominated by neutrinos at the Fermi sur- 
face [35| thus we set, in Eq. 12-13, E v = [i v . A last com- 
ment about the neutrino mean free path within the mixed 
phase Amp: it is simply a weighted sum of the mean free 
paths of the two phases: 1/Xmp = x/A 9 + (1 — X)/X n . 

Let us discuss now our results by starting from the 
case of the deleptonization of an ordinary neutron star, 
which will constitute our reference model. We fix the 
central pressure of the star p c — 83 MeV/fm 3 , and calcu- 
late the equilibrium configuration for the initial uniform 
profiles of S/N and Yl. We stop the integration of the 
structure equations when the pressure drops to a value 
Po corresponding to a baryon density of roughly half of 
the nuclear saturation density [52j . We obtain a baryon 
number of Nb ~ 2.29 x 10 57 with a corresponding grav- 
itational mass of ~ 1.78M Q . We then start the diffusion 
simulation [531 ] ■ Results for the temporal variations of 
the neutrino density, the baryon density and the pres- 
sure profiles are shown in Figs. 3-5. From Fig. 3 one 
can notice that while we are here adopting many sim- 
plifying assumptions for the transport process, the time 
scale of deleptonization is in the right order of magni- 
tude: after t* ~ 7s the central neutrino density is half 
of its initial value. Similar estimates are obtained with 
the simple treatment of the deleptonization and cooling 
time scales presented in [3lj j and in the simulations of 
[36j (on the other hand, in the more recent one [38[, a 
better treatment of the baryon - neutrino interactions 
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FIG. 7: Temporal evolution of the volume fraction of the 
quark phase for a hybrid star of baryonic number Nb = 2.29 x 



10' 



Initially the star has a core of 2SC phase, a thin layer 



of mixed phase and an outer layer of nucleonic matter. After 
~ 4 s a sizable change of the structure of the star occurs 
due to the formation of the normal quark phase triggered 
by deleptonization. The mixed phase has three components: 
nucleonic matter, 2SC and normal quark matter. 



FIG. 8: Temporal variation of the baryon density as a function 
of the enclosed baryonic number for a hybrid star of baryonic 
number Nb = 2.29 x 10 57 . The thick line corresponds to 
the initial configuration t = and the dashed line to the 
configuration at t — 10s. The thin solid lines stand for the 
intermediate configurations with temporal steps of Is. The 
central density slightly decreases during the first seconds but 
after ~ 4s when the normal quark phase appears, it starts to 
increase. 



leads actually to a deleptonization process which contin- 
ues through most of the cooling epoch, tdei ~ 50s). We 
stop our simulation at t = 10s, since for longer times 
neutrino would start to become non-degenerate and our 
assumptions do not hold anymore. 

Instead, the temporal behavior of the density and the 
pressure is different from what was found in [36l.l38| but in 
qualitative agreement with the results of [37| (see Figs. 3 
and 10 of that paper): during the first seconds of the evo- 
lution, while in the outer layers of the star (a/ATg close 
to 1) the density increases due to the lack of neutrino 
pressure, the pressure and density at the center of the 
star slightly decrease. This is actually a consequence of 
the stiffening of the nucleonic EoS when Yl is reduced: 
even if the contribution to the total pressure of neutri- 
nos decreases during deleptonization, nucleonic matter 
becomes more and more isospin asymmetric and, due to 
the large value of the symmetry energy at high density 
within the TM1 model (see [la]), the overall effect is a 
stiffer EoS. This effect deserves more accurate studies 
by implementing a different choice for the initial profiles 
of Yl and S/N taken from supernova simulations and 
by improving the treatment of the neutrino and energy 
transport (energy transport is not included here). Inter- 
estingly, in the recent sophisticated calculation presented 
in Ref. [40( , the central density is almost constant during 
the long term evolution (see Fig . 16a of that paper) at 
variance with the results of [381 ] . Here we want just to 
stress that the symmetry energy could play an important 
role in protoneutron star evolution and, since quark mat- 
ter has usually a smaller symmetry energy than nucleonic 
matter, hybrid stars show a different temporal evolution 
as we will show in the following. 



By looking at the pressure profiles in Fig. 5, we notice 
that the value of the pressure at a = Nb is the same at 
all times as it should be, which represents an important 
check of our calculation. 

Let us discuss now the results obtained for the hy- 
brid EoS presented in Sec. II. We consider an hybrid 
star having the same baryon number as the neutron star 
considered before. The corresponding central pressure 
is p c = 118 MeV/fm 3 and the initial gravitational mass 
AI = 1.78M Q . After calculating the initial configuration 
with constant S/N and Yl, we start the temporal evolu- 
tion simulation. The temporal variation of the neutrino 
density is shown in Fig. 6: notice that the neutrino den- 
sity profile follows the same shape of the baryon density 
profile (Fig. 8), with the typical features of a first order 
phase transition (a mixed phase and two pure phases). 
In the pure 2SC quark phase, the density of neutrinos is 
larger than in the pure nucleonic phase, the central neu- 
trino density is in this case n Vo — 3.44 x 10~ 2 fm~ 3 . This 
is due to the depletion of electrons in the quark phase, 
with respect to the nucleonic phase, and the condition 
of fixed lepton fraction [43l j . The deleptonization time of 
the hybrid star is comparable with the one of the neu- 
tron star: although the mean free path within the quark 
phase is larger than the one in the nucleonic phase at the 
same value of baryon density [44j , neutrinos enter the nu- 
cleonic phase after diffusing inside the quark phase and 
the mixed phase. Since they are more abundant than in 
the case of the neutron star, their mean free path within 
the nucleonic phase is strongly suppressed (see the de- 
pendence of X n on E v in Eq. (12)). The modification of 
the structure and the composition of the star is shown 
in Fig. 7 where the quark volume fraction profile is dis- 



played. This figure represents the main result of this pa- 
per. Starting with an initial configuration with a core of 
2SC phase, a small region of mixed phase and a nucleonic 
matter layer, the star gradually modifies its structure and 
composition during deleptonization, and at some point 
(t ~ 4s in this example), close to the onset of the mixed 
phase, the 2SC pairing pattern is broken and the normal 
quark phase is formed. The mixed phase has now three 
components and occupies a much larger volume fraction 
of the star and continues to grow until a stationary con- 
figuration is reached when almost all neutrinos leak out 
of the star. Quite remarkably, in a short amount of time, 
between 4 and 5 s the structure of the star changes siz- 
ably with possible observable effects. This change of the 
structure of the star is also evident from the baryon den- 
sity profile shown in Fig. 8. The baryon density slightly 
decreases during the first seconds, similarly to the neu- 
tron star case, but after ~ 4s, in coincidence with the 
formation of the normal quark phase in the mixed phase, 
it starts to increase. Initially only a small fraction of 
the star is occupied by the 2SC phase and therefore the 
evolution of the system is dominated by the nucleonic 
phase. With the formation of the normal quark phase, 
very quickly a large fraction of the star is occupied by the 
quark and the mixed phase and the density starts to in- 
crease further favoring the formation of the quark phase. 
This behavior is compatible with the fact that the EoS of 
the mixed phase and the pure quark phase softens as the 
lepton fraction decreases (see Fig. 2.). This qualitative 
difference between neutron stars and stars containing an 
exotic phase was already found in [33] for the case of 
hyperonic matter. Although very interesting, because 
related with the different symmetry energy of the nucle- 
onic phase and the quark phase, the different evolution 
of the density in neutron stars and hybrid stars must be 
investigated in more sophisticated transport models and 
by using more realistic initial conditions before a firm 
conclusion can be drawn. In the same figure, one can 
also notice how the profile of the density qualitatively 
changes during the evolution because of the sizable mod- 
ification of the mixed phase. Unfortunately, we cannot 
provide here signatures of this structural and chemical 
changes in the star within the neutrino signal released; 
such a calculation represents an important extension of 
this work. 

Let us consider one more case with a smaller pro- 
toncutron star total baryon number. Fig. 9 shows the 
temporal evolution of the neutrino density for a neutron 
star having an initial gravitational mass of 1.38M0 and 
Nb = 1-74 x 10 57 . The initial central density is smaller 
than the density of the onset of the mixed phase, the evo- 
lution proceeds as in the case of neutron stars with a slow 
decrease of the central density and after t ~ 6s, when a 
large part of neutrinos already left the star, the normal 
quark phase starts to appear in the mixed phase and 
reaches an equilibrium volume fraction of ~ 0.2, shown 
in the insert of Fig. 9 (this scenario is similar to the 
cases considered in [43] ) . In this case, the 2SC phase 




FIG. 9: Temporal variation of the neutrino density as a func- 
tion of the enclosed baryonic number for a hybrid star of bary- 
onic number Nb = 1-74 x 10 57 . The initial central neutrino 
density, n vo = 1.66 x 10~ 2 fm~ 3 . The thick line corresponds 
to the initial configuration t = and the dashed line to the 
configuration at t = 10s. The thin solid lines stand for the 
intermediate configurations with temporal steps of Is. The 
insert shows the quark volume fraction profiles for the latest 
4 seconds of the evolution. Only the normal quark phase is 
formed in this case at late stages of the deleptonization pro- 
cess. 



cannot be formed because the initial conditions are such 
that quark matter appears when matter is already very 
isospin asymmetric. 

For very large initial total baryon number, a collapse 
is obtained during the evolution because the maximum 
mass is reached. The scenario in that case is quite simi- 
lar to the one obtained for hyperonic matter, kaon con- 
densed matter and normal quark matter |37H39l l42t l43j . 
Due to our simple treatment of the initial conditions and 
the neutrino transport we cannot provide here the evolu- 
tion of configurations close to the maximum mass since 
for those cases the initial conditions arc crucial for the 
stability of the star with respect to collapse. A full calcu- 
lation from the supernova explosion to the late evolution 
of protoneutron stars is needed for very massive stellar 
objects as done in (iol. |4lj|. 



IV. CONCLUSIONS 

We have computed the temporal evolution of pro- 
toneutron stars containing a color superconducting phase 
in the core. We considered here the case of the 2SC 
phase which, between the different color superconducting 
phases, is the most likely to appear in a newly born neu- 
tron star. During the first few seconds after the birth of 
the star, the gradual deleptonization is accompanied by 
a reduction of the electron and proton fractions: matter 
becomes more and more isospin asymmetric. In the 2SC 
quark phase, when the isospin asymmetry reaches a cer- 
tain threshold the pairing cannot take place anymore and 
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a first order phase transition to the normal quark phase 
occurs. Interestingly, the unpairing transition takes place 
close to the onset of the nucleonic matter - quark mat- 
ter mixed phase, where the mismatch between the quark 
chemical potentials is the highest. The structure and 
composition of the star are strongly modified during the 
phase transition and they reach an equilibrium configura- 
tion at late times when most of the neutrinos are released. 
This is the typical evolution of a star with initial grav- 
itational mass larger than ~ 1.6M . For less massive 
hybrid stars, quark matter appears only at later times 
and the isospin asymmetry is then so large that the 2SC 
phase cannot be formed and only a normal quark phase 
can appear in the core of the star. 

The calculations we have presented here are based on 
many simplifying assumptions for the neutrino interac- 
tions and for the treatment of the neutrino transport. 
We have nevertheless provided a few examples of how the 
evolution of a protoneutron star can be rich and com- 
plex when color superconducting phases are taken into 
account. More realistic calculations are of course needed 
and also the three flavor CFL phase should be included 
in the equation of state. Work along this line is already 
in progress. 

Finally we want to remark why the scenario we have 
discussed here is potentially very interesting for phe- 
nomenology: the neutrino signal, could be actually mod- 
ified by the happenings of the core of the star as shown 
in p5l - l47| for the scenario of a phase transition to quark 
matter in the early post-bounce stage of a supernova. 
If sizable modifications of the structure of the star oc- 
curs one could indeed expect to have signatures on the 



spectrum and the luminosity of the emitted neutrinos: 
for instance, the average energy of the emitted neutrinos 
strongly depends on the compactness of the star [40| . 
Also the effect of metastability and nucleation connected 
with a first order phase transition [48| could have an im- 
print on the temporal structure of the signal. The present 
neutrino detectors, as SuperK, will be able to catch thou- 
sands of neutrinos for a galactic supernova thus allowing 
a detailed temporal and spectral investigation of the sig- 
nal m. 

Another very exciting implication of our scenario con- 
cerns gamma-ray-bursts: there are now many hints on 
the possibility that at least some gamma-ray-bursts are 
produced by strongly magnetized neutron stars. The 
magnetar model of gamma-ray-bursts (49| as compared 
to the collapsar model could explain better the long term 
activity, up to 10 6 s after the prompt emission, seen in 
some light curves [5(3, [HJ • This implies that the prompt 
emission is connected with the first stages of the evolu- 
tion of a protoneutron star and the specific features of 
the neutrino signal could have therefore an impact on 
the gamma-ray-bursts. The complex temporal structure 
of some gamma-ray-bursts could have its explanation in 
the complex phenomena which can occur in a protoneu- 
tron star if quark phases are formed. 
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